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Abstract. The paper is concerned with the finite element solution of the Poisson equation with ho- 
mogeneous Dirichlet boundary condition in a three-dimensional domain. Anisotropic, graded meshes 
from a former paper are reused for dealing with the singular behaviour of the solution in the vicinity 
of the non-smooth parts of the boundary. The discretization error is analyzed for the piecewise linear 
approximation in the i/ 1 (0)- and L 2 (fi)-norms by using a new quasi-interpolation operator. This new 
interpolant is introduced in order to prove the estimates for L 2 (f2)-data in the differential equation which 
is not possible for the standard nodal interpolant. These new estimates allow for the extension of certain 
error estimates for optimal control problems with elliptic partial differential equation and for a simpler 
proof of the discrete compactness property for edge elements of any order on this kind of finite element 
meshes. 
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1. Introduction 

We consider the homogeneous Dirichlet problem for the Laplace equation, 

- Au = f inn, u = on 8Q, (1.1) 

where £1 is a polyhedral domain. Note that we could consider a more general elliptic 
equation of second order. But by a linear change of the independent variables the main 
part of the differential operator could be transformed to the Laplace operator in another 
polyhedral domain such that it is sufficient to consider the Laplace operator here. 
The aim of the paper is to prove the discretization error estimate 

\\u - Ufc||ffi(n) < Ch\\f\\ L 2 {n) (1.2) 
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for the finite element solution Uh E Vh which is constructed by using piecewise linear 
and continuous functions on a family of appropriate finite element meshes Th- Note that 
we assume here not more than / E L 2 (Q) such that the L 2 -error estimate 

\\u - u h \\ L 2 {n) < C7i 2 ||/|| L 2 (Q) (1.3) 

follows by the Aubin-Nitsche method immediately. The generic constant C may have 
different values on each occurrence. 

If the solution of the boundary value problem (1.1) was in H 2 (Q) then the finite el- 
ement meshes could be chosen quasi-uniform, and the error estimates (1.2) and (1.3) 
would be standard. However, if the domain Q is non-convex, the solution will in general 
contain vertex and edge singularities, that means u H 2 (ft). In this case the conver- 
gence order is reduced in comparison with (1.2) and (1.3) when quasi-uniform meshes 
are used. As a remedy, we focus here on a priori anisotropic mesh grading techniques 
as they were investigated by Apel and Nicaise in [4]. In comparison with isotropic local 
mesh refinement, the use of anisotropic elements avoids an unnecessary refinement along 
the edges. 

The estimate (1.2) is in general proven by using the Cea lemma (or the best approxi- 
mation property of the finite element method), 

\\u-u h \\ H i {n) <C inf l|tt-«/»||fri(n)> (1-4) 

and by proving an interpolation error estimate as an upper bound for the right-hand 
side of (1.4). The particular difficulty is that when the Lagrange interpolant is used 
together with anisotropic mesh grading, then the local interpolation error estimate 

\u - Ihu\ w i,p( T } < hr\u\ W 2,p(T) (1.5) 

does not hold for p = 2 but only for p > 2, see [2]. Hence the classical proof of a finite 
element error estimate via 

II" - u h \\ H i( Q) < C\\u - I h u\\m(n) < C I ^2 h T \u\ 2 H2(T) 

\TeT h 

does not work. This problem was overcome by Apel and Nicaise, [4], by using (1.5) and 
related estimates in weighted spaces, as well as the Holder inequality for the prize that 
/ G L p (£l) with p > 2 has to be assumed in problem (1.1). Hence estimate (1.2) cannot 
be proved in this way. 

For prismatic domains and tensor product type meshes the problem was overcome 
in [1, 6] by proving local estimates for a certain quasi-interpolation operator. This 
work cannot be easily extended to general polyhedral domains since the orthogonality 
of certain edges of the elements was used there. The aim of the current paper is to 
construct a quasi-interpolation operator such that the error estimate 

\\u - D h u\\ H i(n) < Ch\\f\\ L 2 {n) (1.6) 



/ 
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can be proved for the anisotropic meshes introduced in [4]. 

Quasi- interpolants were introduced by Clement [14]. The idea is to replace nodal 
values by certain averaged values such that non-smooth functions can be interpolated. 
This original idea has been modified by many authors since then. The contribution by 
Scott and Zhang [30] was most influential to our work. 

The plan of the paper is as follows. In Section 2 we introduce notation, recall regularity 
results for the solution u of (1.1) and describe the finite element discretization. The 
main results are proved in Section 3. The paper continues with numerical results in 
Section 4 and ends with two sections where we describe applications which motivated 
us to improve the approximation result from \\u — UhW^tn) — C^||/||iP(Q)j P > 2, to 
|| u — Uh\\ H ifQ\ < Ch\\f\\ L 2(Qy The first one is a discretization of a distributed optimal 
control problem with (1.1) as the state equation. The second application consists in a 
simpler proof of the discrete compactness property for edge elements of any order on 
this kind of finite element meshes. 

We finish this introduction by commenting on related work. The idea to treat sin- 
gularities due to a non-smooth boundary by using graded finite element meshes is old. 
The two-dimensional case was investigated by Oganesyan and Rukhovets [25], Babuska 
[9], Raugel [27], and Schatz and Wahlbin [28]. In three dimensions we can distinguish 
isotropic mesh grading, see the papers by Apel and Heinrich [3] and Apel, Sandig, and 
Whiteman [5], and anisotropic mesh grading, see the already mentioned papers [2, 1, 6] 
for the special case of prismatic domains, and [4] for general polyhedral domains. This 
work has been extended by Bacu^a, Nistor, and Zikatanov [12] to higher order finite 
element approximations where naturally higher regularity of the right-hand side / has 
to be assumed. Boundary element methods with anisotropic, graded meshes have been 
considered by von Petersdorff and Stephan [26]. The main alternative to mesh grading 
is augmenting the finite element space with singular functions, see for example Strang 
and Fix [32], Blum and Dobrowolski [11], or Assous, Ciarlet Jr., and Segre [8] for various 
variants. It works well in two dimensions where the coefficient in front of the singular 
function is constant. In the case of edge singularities this coefficient is a function which 
can be approximated, see Beagles and Whiteman [10], or it can be treated by Fourier 
analysis, see Lubuma and Nicaise [22]. 

2. Notation, regularity, discretization 

It is well known that the solution of the boundary value problem (1.1) contains edge and 
vertex singularities which are characterized by singular exponents. For each edge e, the 
corresponding leading (smallest) singular exponent A e is simply defined by A e = ir/uj e 
where uj e is the interior dihedral angle at the edge e. For vertices v of D,, the leading 
singular exponent A^ > has to be computed via the eigenvalue problem of the Laplace- 
Beltrami operator on the intersection of f2 and the unit sphere centered at v. Note that 
A e > | and A t , > 0. A vertex v or an edge e will be called singular if A„ < i or A e < 1, 
respectively. We exclude the case that ^ is & singular exponent of any vertex. For a 
detailed discussion of edge and vertex singularities we refer to [16, Sections 2.5 and 2.6]. 
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As in [4] we subdivide the domain f2 into a finite number of disjoint tetrahedral 
subdomains, subsequently called macro -elements, 



n 



We assume that each Ae contains at most one singular edge and at most one singular 
vertex. In the case that Ae contains both a singular edge and a singular vertex, that 
vertex is contained in that edge. Note that the edges of Ae are considered to have O(l) 
length. For l\ ^ £2, the closures of the macroelements and Ag 2 may be disjoint or 
they intersect defining a coupling face, or a coupling edge, or a coupling node. Denote by 
T c , £ c and M c the sets of coupling faces, edges and nodes, respectively. 

For the description of the regularity of the solution u of (1.1), we set A^ = A^ if 
the macro-element Ae contains the singular vertex v of 0. If Ae does not contain any 



singular vertex we set \i = +00. Moreover, we set A e = A e if Ae contains the singular 

(£) 

edge e of f2, otherwise we set Ae = +00. Furthermore, we define in each macro-element 



(0 



Ae a Cartesian coordinate system x 



(f) 



x 



1 ) x 2 ' x 3 



such that the singular vertex, if 



existing, is located in the origin, and the singular edge, if existing, is contained in the 



axis. We also introduce by 



0W{xW) : 



( { x^? + { xf?) 1/2 , 
):= {^f + ( X f? + { X ff) ll \ 

RW(x®y 



(£) 

the distance to the x% -axis, the distance to the origin, the angular distance from the 
respectively. 

For k G N and j3, 5 G K we define the weighted Sobolev space 

y";I(^) ■= e V'{Ai) : NlyM (A< ) < 00 } 



where 



^2 f RP- k+ \ a \0 s ~ k+ \ a \D a i 

a\<k^ Ae 

[ R P 5 D a t 
_i_,<VA» 



\a\=k Uiyt 



Here, we have used the standard multi-index notation to describe partial derivatives, 
and we have omitted the index (£) in R and 6 for simplicity. 
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Theorem 2.1 (Regularity). [4, Theorem 2.10] The weak solution u of the boundary 
value problem (1.1) admits the decomposition 



u = u r + u s 



in At, £ = 1, . . . , L, where u r S H 2 



(Ae) and 





for any (3,5 > satisfying f3 > | — Av end 5 > 1 — \ e . 
Following [4] we consider a triangulation Th of O, 



n= U r ' 



made up of tetrahedra which match the initial partition: if T n 7^ then T C A^. 
Four cases are considered: 

1. If A^ does neither contain a singular edge nor a singular vertex then Th\ A e is 
assumed to be isotropic and quasi- uniform with element size h, see Figure 1, top 



2. If A^ contains a singular vertex but no singular edges then Th\A e is isotropic and has 
a singular vertex refinement, i.e., the mesh is graded towards the singular vertex 
with a grading parameter vi G (0,1]. This can be achieved by using a coordinate 
transformation of the vertices from Case 1, see Figure 1, top right. 

3. If Ai contains a singular edge but no singular vertices then Th\h t is anisotropically 
graded towards the singular edge. The grading parameter is im G (0,1]. To 
this end, we introduce a family Vi of planes transversal to the singular edge and 
containing the opposite one. These planes split the macro element into strips and 
contain all nodes. In the planes the position of the nodes is achieved by applying 
a coordinate transformation to a uniform triangulation, see Figure 1, bottom left. 

4. If Ai contains both a singular vertex and a singular edge then Th\\ e is graded 
towards the singular edge with grading parameter [ii 6 (0, 1] and towards the 
singular vertex with grading parameter G (0,1]. The mesh is topologically 
equivalent to the mesh of Case 3 but the planes of Vi do not divide the singular 
edge equidistantly but with a grading towards the singular vertex. 

We point out that anisotropic elements can appear only in Cases 3 and 4, for which 
Th contains needle elements near the singular edge and flat elements near the opposite 
one, see Figure 1. We further observe that if Ag is of type 3 or 4, the elements in Th\\i 
do not intersect any plane of Vi- 



left. 
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Figure 1: Macroelements of types 1, 2, 3 and 4 



For each element T we introduce its lengths hi t T, h,2,T, ^3,T and as follows. Let hx 
be the diameter of T. If T C with of type 1 or 2, then /iit = /i2T = ^3,T = /it- 
If T C A^ with A^ of type 3 or 4 then h^ t T is the length of the edge e^p of T parallel to 
the singular edge, and = h^T = + l e 2,T|) where e\g and e2,T are the edges 

of T intersecting e% t and each one of them is contained in some plane of Vi- 

By classical regularity theory, the solution u of the boundary value problem (1.1) is 
continuous, see e.g. [16, page page 79], such that the Lagrange interpolant uj with 
respect to the subdivision {A^} is well defined. We consider the decomposition 

u = uj + ur. (2-1) 

It follows that the restriction ur\a ( has the same smoothness properties as u, see The- 
orem 2.1. Furthermore, ur vanishes in coupling nodes and on singular edges. We 
construct now an interpolant D^ur S Vh which also vanishes on these nodes such that 
ui + D^ur £ Vh can be used to estimate the discretization error via (1.4). 

To this end, let M, A/i n , M c and M s be the set of all nodes of 7~h, the set of all the 
interior nodes, the set of coupling nodes, and the set of nodes which belong to some 
singular edge, respectively. The terminal points of the singular edges are included in 
M s . The piecewise linear nodal basis on Tk is denoted by {4> n }neJV- We associate (as 
specified below) with each n G M \ (M c U M s ) an edge a n with n as an endpoint. Note 
that u\ an 6 L 2 (a n ) since u £ H S (Q) with s > 1. Hence the operator Dh with 

D h u= Yl (IWi)(n) • 0„(a;), (2.2) 
neAf\(M c uM s ) 
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Figure 2: Illustration of the edges a. 



is well denned when : L 2 (a) — > V\(o) is the L 2 (cr)-projection operator onto the space 
of polynomials of degree less than or equal to one. Note that D^u vanishes on coupling 
nodes and on singular edges by construction. In order to impose the boundary conditions 
and to be able to prove interpolation error estimates we need to select the edges a n in 
an appropriate way, compare the illustration in Figure 2. First, we demand that 

• for each node n G N \ (A/" c U A/" s ), n and a n belong to the same macroelement. 
This requires in particular the following restrictions. 

— If n lays on a boundary or coupling face, then a n is contained in that face. 

— If n lays on a coupling edge, then a n is contained in that coupling edge. 

Note that these requirements made the treatment of the coupling nodes via the in- 
terpolation on the initial ui necessary. Note further that this construction leads to a 
preservation of the homogeneous Dirichlet boundary condition. 

In order to prove the stability of in the anisotropic refinement regions we also 
require: 

• If n is a vertex of a tetrahedron contained in a macroelement Ag of types 3 or 4, 
then o n is an edge contained on some plane of Vt. 

• If n\ and belong to a macroelement Ai of types 3 or 4 and have the same 
orthogonal projection onto the x\ x\ -plane, then the same holds for a ni and a n2 . 

In order to estimate the interpolation error we need to define for each T £ Th a set 
St which should satisfy the following assumptions. 

• The set St is a union of elements of Th (plus some faces) and in particular T St- 

• The set St is an open connected domain, and as small as possible. 

• We have a n C St for all nodes n of T. 

• If T C A e , then S T C A e . 

• If T C Ai with Ai of type 3 of 4, then St is a prism where the top and bottom 
faces are contained in two planes of (and so they are not parallel) and the other 
faces are parallel to the singular edge. 
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The following properties follow from the definitions of the edges a n and the sets St- 

1. Let T be contained in a macroelement of type 3 or 4. If T intersects two planes 
pi and p2 of Vii then St intersects exactly the same planes p\ and p2- 

2. If the node n, n Af c U A/" s , belongs to a coupling face, that means that there exist 
tetrahedra T\ C A^ and T2 C A^ 2 with i\ ^ £2 and n E R T2, then S^i R St 2 = 
but cr n C Stj R $t 2 • 

3. If T is an isotropic element then all the elements in St are also isotropic and of 
size of the same order. 

The second point is essential for our proof of the approximation properties. It was the 
target for which we made the construction as it is. 

3. Error estimates 

The aim of this section is to derive error estimates for our discretization. They are 
based on local interpolation error estimates for our interpolant D^. For proving these 
estimates we have to distinguish several cases, see also Figure 3 for an illustration: 

1. T is an isotropic element without coupling node, u has full regularity, 

2. T is an isotropic element with coupling node, u has full regularity, 

3. T is an isotropic element with coupling node, u has reduced regularity, 

4. T is an anisotropic flat element without coupling node, u has full regularity, 

5. T is an anisotropic flat element with coupling node, u has full regularity, 

6. T is an anisotropic needle element without node on the singular edge, u has full 
regularity, 

7. T is an anisotropic needle element with node on the singular edge, u has reduced 
regularity. 

In Lemma 3.1 we present the general approach for the proof of the local interpolation 
error estimate by considering isotropic elements with and without coupling nodes (cases 
1 and 2). We proceed with Lemmas 3.2 where we introduce for isotropic elements how to 
cope with the weighted norms in the case of reduced regularity (case 3). The interpolated 
function is only from a weighted Sobolev space but we will see that this even simplifies 
some parts of the proof. 

For anisotropic elements the use of an inverse inequality (as was done in the pre- 
vious lemmas) has to be avoided; instead we use the structure of the meshes in the 
macroelements of types 3 and 4. We start with a stability estimate of d^D^u which 
allows immediately the treatment of anisotropic flat elements (cases 4 and 5) in Lemma 
3.4. Then we prove stability estimates for the remaining derivatives and continue with 
the interpolation error estimates for needle elements. Lemma 3.7 is devoted to case 6, 
and Lemma 3.9 to case 7. 
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Figure 3: Illustration of the cases that have to considered for the interpolation error 
estimates 



All these local estimates can then be combined to prove the global interpolation error 
estimate, see Theorem 3.10, and the finite element error estimates, see Corollary 3.12. 

Lemma 3.1 (isotropic element, full regularity). If T is an isotropic element then the 
local interpolation error estimate 

\u - D h u\ H i {T) < Ch T \u\ H 2 (ST) (3.1) 

holds provided that u £ H 2 (St) and u(n) = for all n G M c . 

Proof. Following the explanations in [30, page 486] and [1, page 1156], an explicit repre- 
sentation of D^u from (2.2) can be given by introducing the unique function ip n £ Vh\a n 
with J a tp n (f>j = 5 n j for all j E M such that 

(n CTn n)(n)= ! u^ n (3.2) 

J On 

and 



D h u\ T = XM / u ^ n ) ' ^ n 



(3.3) 



neA/r 

where we denote by Mt the set of nodes of T without the coupling nodes. Note that 



t] n\\L°°{a n ) = C\a n \ , (3.4) 
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compare [1, page 1157]. (By some calculation one can even specify that C = 4.) With 
(3.3), the direct computation 

\^n\m(T) < Ch^}\T\^ 2 , (3.5) 

the trace theorem 

IMUi(0 < C\a n \\S T \^ ll2 {\\u\\ L 2 {ST) + h T \u\ H i (ST) + h 2 T \u\ H 2 [ST) ), (3.6) 
and \St\ < C\T\ we obtain 

\DhU\ H i(T)<C ^ ||^||L 1 ( ( r n )IIV'n||L°°( ( T n )l0n|// 1 (T) 
neA/r 

< Ch^iWuW^^) + h T \u\ m{ST) + h 2 T \u\ H 2 {ST) ). (3.7) 

If Mt does not contain a node n G M c we find that D^w = to for all w G Vi such that 
we get by using the triangle inequality and the stability estimate (3.7) 

\u - D h u\ H i( T ) = \(u-w)- D h (u - io)|fl-i(T) G Pi 
<\u- w\ H i( T ) + - w)\h^(T) 

< C (h^Wu - w\\ L 2 (St) + \u- w\ h i {St) + h T \u\ H 2 {ST) ) . 

We use now a Deny-Lions type argument (see e.g. [15]) and conclude estimate (3.1). 

In the case when Mt contains a node n G N c we do not have the property that 
DhW = w for all w G V\ but we can use that u(n) = 0. Let a n be an edge contained 
in T having n as an endpoint, and let <p n be the Lagrange basis function associated 
with n. (Note that we deal here with nodes n which are not used in the definition of 
Dfr. Therefore we can assume that a n is local in A^.) Consequently, we have with the 
previous argument that 

\u - (D h u + (n ff „u)(n)0 n )| Hl(T j < Ch T \u\ H 2(s T y (3.8) 

Let Itu be the linear Lagrange interpolation of u on T. Since lTu\a n is linear, we have 
have (Tl an lTu){n) = 0. Prom this fact and using (3.2)-(3.6) as in the derivation of (3.7) 
(here with the specific T instead of St since a n C T), we have 

|(n <Tn M)(n)0 n | H i (T) = KU^iu - I T u))(n)4> n \ H i {T) 

< Ch T l (ju — Itu\ L 2( T ) + Ht\u — Itu\hi(T) + ^t\ u \h 2 (T)) 

< Ch T \u\ H 2( T ) 

where we used standard estimates for the Lagrange interpolant in the last step. With 
(3.8) and the triangle inequality we conclude estimate (3.1) also in this case. □ 

Lemma 3.2 (isotropic element, reduced regularity). If T is an isotropic element then 
the local interpolation error estimate 

\u - D h u\ H i( T ) < Ch T ~P\\u\\ v 2,2^ ST) (3.9) 



2 2 

holds provided that u G Vg ' (St), ft G [0, 1) 
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Proof. We start as in the proof of Lemma 3.1 but use the sharper trace theorem 

IMUi(<r„) < CWnWSrl^iWuWLi^ + h T \u\ w iA {ST) + h 2 T \u\ W 2,i {ST) ). 
With (3.3), (3.4), (3.5), and \S T \ < C\T\ we obtain 

\D h u\ H l(T) < C \\u\\L^(a n )\\^n\\L^(u n )\^n\m{T) 
raGA/r 

< CI-S'Tr^^lhlk 1 ^) + Hw^(St) + h T\u\w^{s T )) 

< C(^ 1 [|«|| £ a( S i r ) + Mhi(St) + \ S T\' 1/2 h T \u\ W 2,i { s T) ) 

and hence via the triangle inequality 

\u - D h u\ H i (T) < C(^ 1 [|«|| £ a( S i r ) + Mhust) + \ST\~ 1/2 h T \u\ W 2,i(s T) ). (3.10) 
For the first two terms we just use that R < hx, hence 1 < hxR~ l , to get 



|«|tfi(S T ) < 4 



1-/3 



n ,,i,2 



^-i,o(St)- 



To estimate the third term we use the Cauchy-Schwarz inequality and again R < hx,to 
obtain for lad = 2 



< C\S T \ 1/2 h T -\v\ t 



where \\R ^\\l 2 {s t ) — C|5t| 1/ ' 2 /i t /3 is obtained by executing the integration and using 
that /3 < |. All these estimates imply estimate (3.9). □ 

In order to prove interpolation error estimates for the anisotropic elements we derive 
stability estimates for where we avoid the use of the inverse inequality. Let x\,x<i 
and X3 be a Cartesian coordinate system with the ^-direction parallel to the singular 
edge of A. We will estimate separately the L 2 -norm of the derivatives of D^u. 

Let T be an anisotropic element with the characteristic lengths = h^T an d h^T- 
We will not use that h^T > ^j,T> j = 1, 2, in the next lemma in order to use this estimate 
both for the needle and the flat elements. 

Lemma 3.3 (Stability in direction of the singular edge). For any anisotropic element 
T the estimate 

\\d 3 D h u\\ L2(T) < c\St\- 1/2 Yl mD a d 3U \\ms T) 

\a\<l 

holds provided that 8311 S W 1 ' 1 (St)- 
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Proof. We observe that T has an edge ex parallel to the singular edge, and so, parallel to 
the X3-axis. Since D^u is linear on T, we have d 3 Dfru\T = d 3 Dhu\ eT . If ex is contained 
on the singular edge, then d 3 Dhu\T = since Dhu\ eT = u\ eT = and we are done. Now, 
consider the case that ex is not contained in a singular edge and denote its endpoints 
by m and ri2 such that d 3 cj) ni \ T = -/i 3 y and d 3 <fi n2 \ T = h 3T . Then we have 

d 3 D h u = 

We observe now that by our assumptions cr ni and a n2 have the same projection o~t 
into the xi^-plane and hence form two opposite edges of a plane quadrilateral which is 
parallel to the X3-axis and which we will denote by Ft- We note further that ip ni and ip n2 
can be considered as the same function Tpx defined on cry and ||V't||l° o (<tt) = C\o-t\~ 1 - 
With this insight we obtain 



d 3 D h u\ = h 3 1 


/ nV>n 2 - / WVu 


= h 3 l 


/ d 3 uip T 











<Ch 3 \a T \ \\d 3 u\\ L i (FT) <C\F T \ \\d 3 u\\ L i {FT y 



We integrate this estimate over T, apply the standard trace theorem 
\\v\\ lHFt) ^CIFtWStI- 1 h T \\D a v\\ L x {ST) 

\a\<l 

and obtain the desired estimate. □ 

We are now prepared to estimate the interpolation error for the flat elements occurring 
far away from the singular edge in cases 3 and 4. 

Lemma 3.4 (anisotropic flat element, full regularity). IfT is an anisotropic flat element 
(h 3t T < h±T = ^-2,T ) then the local interpolation error estimate 

\u - D h u\ H i( T ) < Ch T \u\ H 2( ST) (3.11) 

holds provided that u £ H 2 (St)- (Remember that h? = diam(T). / ) 

Proof. The proof for d 3 (u — D^u) can be done on the basis of Lemma 3.3. Assume for 
the moment that the element T does not contain a coupling node. Similar to the proof 
of Lemma 3.1 we obtain for any w & V\ 

||d 3 (tt - DhV^Wtf^ = \\d 3 (u - w) - d 3 D h (u - w)\\ L 2^ 

< C\\d 3 (u - w)\\ lHSt) + C ]T h T \\D a d 3 u\\ L i {ST y 

M=i 
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We choose now w G V\ such that the constant d 3 w satisfies J s d$(u — w) = and such 
that we can conclude by using the Poincare-Friedrichs inequality (or again a Deny-Lions 
type argument) 

\\d 3 (u-w)\\ LHST) <C^ h$\\D a d 3 u\\ LHST) 

H=i 

and hence 

\\d 3 (u - D h u)\\ L2{T) <C h T\\ Da dM\^(s T ) < Ch T \u\ H 2 ( s T) . 

\a\=l 

Note that the polynomial w can be chosen such that it vanishes in three nodes of T. It 
is completely described by choosing the appropriate value at one endpoint of the edge 
of T which is parallel to the ^3-axis. Since a possible coupling node is not an endpoint 
of this edge, the argument above can also be used in the case of coupling nodes. 

For the other directions we can proceed as in the proof of Lemma 3.1. In the case 
of coupling nodes the interpolation error estimate \u — Itu-Ih 1 ^) — C^rM/p^) is used 
there which does not hold for anisotropic elements. However, the estimate \\di(u — 
Itu)\\l2(t) < CIit\u\ H 2( T }, i = 1,2, does hold, see for example [2]. □ 

It remains to prove interpolation error estimates for needle elements such that we will 
assume hi x = hi.r < C^3,r for the next lemmas. 

Lemma 3.5 (Stability in direction perpendicular to singular edge, anisotropic needle 
element away from singular edge). Assume that the element T does not contain a node 
n G M s and that h\ t T = ^2,T < Ch% y T- Then for i = 1,2 we have 

\\diD h u\\ L 2 {T) < C {\u\ h i {St) + h^ T \d 3 u\ H i {ST) ) (3.12) 

provided that u S H (St) o,nd d 3 u G H 1 (St)- 

Proof. For each node n S Mr we denote by F n ^r the top or bottom face of the prismatic 
domain Sr such that n £ F Hi r- Observe that we have a n C F n: r C St for all n G Mr- 
Observe further that F Ui t is isotropic with diameter of order h\ t T and recall the standard 
trace inequality 

lb|| L i ((Tn) < C^nWF^- 1 (\\v\\ LHFnT) + h^ T \v\ w i,i {FnT) ^ (3.13) 
for all v G W 1 ' 1 (F n ^r)- We need also the trace inequality 

[|«[Ui(F n , T ) < C^tWStI 1 (\\v\\ lHSt ) + ^3,T||^|| L i (5T )) (3.14) 

which can be proved by using Lemma A.l from page 23 and the facts that St is a union 
of prisms, and F n p< is a face of St- 

Let st be one of the short edges of T and denote its endpoints by n 1 and n 2 . We use the 
same notation sr for the direction of this edge in order to denote by d ST v = Vv • st/\st\ 
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the directional derivative. In the following we first estimate \\d ST Dhu\\i2^ T y After that, 
the desired estimates (3.12) easily follow as we will show. 

Notice that if n € ATt \ {n 1 ,?! 2 } we have d ST (j) n = 0, and if n S {ra 1 ,?! 2 } then 
l|^ST0n|U°°(T) = I s t| _1 < Ch^ T . For all w € Pq(St) we have (and here we use that the 
element does not contain a node n 6 Af c UA/" S ) 



\d ST D h u\\ L 2( T) = \\d ST D h (u- 



I ft 



WvMlL 2 (T) 



^Ch^T] 1 / 2 \\ u - w \\Li(ar n )\\'>Pn\\L°°(a n ) 
< C^ITI 1 ^ J2 lan^Wu-wW^^y (3.15) 
From the trace inequality (3.13) we have for each n S Mt D st 

HI^K) - Cknll^rl" 1 - w|| L i(F nT ) +hi )T \u\ w i,i {Fn T) 



\\u 

Since the definition of F n ^ implies F n i = F n i =: Ft, we have 

II" - w \\^(a n ) < C^nWFr^ 1 [\\u - w\\ l i {Ft) + h 1>T \u\ w i,i( FT) ) . 

Now we choose uu as the average of u on F F and use a Poincare type inequality on Fx 
to get 

Therefore we arrive at 

\\d ST D h u\\ L 2 {T) < C\T\ l l 2 \F T \~ 1 \u\ w \,i(p T - ) 

< CITI^St! -1 {\u\w^{s t ) + h 3,T\d3u\ w i,i(s T )) 

< C\S T \~ 1/2 (\u\ w i,i {St) + h 3tT \d 3 u\ w i,i(s T) ) (3.16) 

< C (|u| h i (St ) + h 3tT \d 3 u\ H i {ST) ) 

where we used again the trace inequality (3.14). 

Now, let si t T and S2,T be two different short edges (edge vectors) of T such that the 
determinant of the matrix made up of , 1,T , , , 2,T i and e.q as columns is greater than a 

1 |«1,t|' |«2,t| ta 

constant depending only the maximum angle of T. Note that this is possible due to the 
maximal angle condition, see [18]. Then, if the canonical vector ej, i = 1, 2, is expressed 
as 

Sl.T $2,T 
ei = Cii, r + C2,i, r + C3,je3, 

l s i,r| |S2,t| 

it follows that ci^, C2,i and C3,i are bounded by above by a constant depending only on 
the maximum angle condition. Since 

di = c lti d SlT + c 2 ,id S2T + c 3ii 5 3 

we obtain (3.12) from (3.16) with st = sit arid st = S2T1 Lemma 3.3, and recalling 
that h\ t T = ^2,T < C^3,T- D 
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Lemma 3.6 (Stability in direction perpendicular to singular edge, anisotropic needle 
element at the singular edge). Assume that the element T contains at least one node 
n G M s and that H\t = ^2,T < C^3,T- Then we have for i = 1, 2 



\diD h u\\ L 2 {T) 



< C\S 7 



-1/2 



u 



W^{S T ) + T^II<93«IUi(S T ) + ^ ^tI^^I^M^) I ( 3 ' 17 ) 

|a|=l 



provided that u G 1U 2,1 (,St)- 

Proof. For each node n G A4 of T we select one short edge a n with an endpoint at n 
and contained in the same macroelement as T such that we can apply Lemma 3.5. We 
have for i = 1,2 



L2( T) 



9i ^2 (n CTn n)(n)^ n 

\ neAfsriT 

< C\S T \~ 1/2 (\u\ w i,i {St) + /i 8 ,T|^«| w n,i(5 T )) . 

Now we deal with ||3i[(n o - n w)(n)(/) n ] 11^2^ which is first estimated by 

II^Kn^uXnJcOlram <CW 11 -' 1-11 



(3.18) 



(3.19) 



|L 2 (T) ^ Ml^nllZPCT)! "™! IFIlL^o-n) 
for each n G A/" s fl T. 

Let n £ M s nT and be i^r be the face of St having a n as an edge and another 
edge on the singular edge. Let P Ht T be the greatest parallelogram contained in F Uj t and 
having a n as an edge. So, P n ,T is parallel to the x 3 -axis, and its area is comparable with 
the area of F n) T since opposite edges of the trapezoid F n ,T have equivalent length. Using 
a trace inequality we have 

- C, | cr ™ll i? n,Tr 1 (||n|| L i(p n T ) + /i 3)T ||«93w|| L l(p niT )). 

But, since u = on the edge of P n ,T contained on the singular edge we can use the 
Poincare inequality to obtain 

< C\a n \\F n ^ T \- l {\cj n \\\d an u\\ L i {PnT) + h 3 , T \\d 3 u\\ L i iPnT) ). (3.20) 
From Lemma A. 2 we have for all v G W 1 ' 1 (St) 

\\ v \\LHP ntT ) 

< CIF^tWStI' 1 (IMIlI(St) + \ s hT\\\ d si,T V \\LHS T ) + \s2,T\\\ d S2,T V \\Li(S T )) ■ ( 3 ' 21 ) 

Using twice (3.21) we obtain from (3.20) 

< Clan^lSrl' 1 (\\da n u\\ L i(s T ) + \ s i,T\\\ d si,T<rn u \\Li(s T ) + |s 2j t|||9 S2T<Tw «|| l i (5t) ) + 
+ C\a n \\S T \~ l h^ T (||<9 3 u|| l i (5t) + |si, T | ||5 S1iT3 m|| l i (St) + |s 2 ,t| ||5 S2t3 u|| l i (5t) ) . 

(3.22) 
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With the estimates 

\\d*n u \\Li(S T ) ^ klw 1 . 1 ^)' 
\\^s i> Tcrn u \\L 1 {S T ) - \ u \w^{S T ), * =1 > 2 , 
H$>i,r3w||zi(S T ) - I^Hw^OSr)' i = 1 > 2 , 

the inequality 

||^0„|| L 2 (T) < C^lf/ 2 , 
and |<r n | ~ /ijT (i = 1,2) we obtain from (3.19) 

||<9i[(n CTn u)(n)0 n ]|| L 2 (T) < C\S T \' 1/2 (\u\ w i,i {St ) + {hx >T + h 2 ,T)\u\ W 2,i { s T) ) 

+ C\S T \- 1 / 2 ^\\d 3 u\\ LHST) . (3.23) 

Finally, taking into account that, since hi T = h 2 ,T < C^3,r, we have 

(h hT + h 2 ,T)\u\ W 2,i(s T ) + h 3 ^\\d3u\\ w i,i {ST) <C h%\D a u\ w i,i {ST ), 

l«l=i 

inequality (3.17) follows from (3.18) and (3.23). □ 

We are now prepared to estimate the interpolation error for needle elements. 

Lemma 3.7 (anisotropic needle element, full regularity). IfT is an anisotropic element 
with hi T = h 2) T < C^3,T then the local interpolation error estimates 

\u - D h u\ m{T) < C h%\D a u\ H i {ST) (3.24) 
l«l=i 

hold provided that u G H 2 (St)- 

Remark 3.8. The estimate (3.24) does not hold for the Lagrange interpolant, see [2]. 

Proof. (Lemma 3.7) Since the needle elements with full regularity do not contain a 
coupling node we can apply both Lemmas 3.3 and 3.5. That means we have shown that 

\D h u\ H i {T) < C h T\\ Da dM\tf{ST) + C {\ u \h\St) + h 3,T\ d 3u\m(S T )) 

H<i 

<C Y h^\D a u\ m{ST y 
\a\<l 

We exploit now that D^w = w for all w G Vi- Consequently, we get 
| it - D h u\ H i(r) = \(u-w)- D h (u - w)\ H i( T ) Vu> G V\ 
<\U- w\ H l(T) + \D h (u - w)\ H l(r) 

<CY h^\D a {u-w)\ m{ST) . 
H<i 

We use now again a Deny-Lions type argument where the form of Lemma 1 in [1] best 
suits our needs, and conclude the desired estimate. □ 
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Lemma 3.9 (anisotropic needle element, reduced regularity). Let T be an anisotropic 
element with hi^ = ^2,r 5; C^3,T and let St have zero distance to the singular edge. 
Then the local interpolation error estimate 

2 

\u - D h u\ H i(T) < Ch\p \\diu\\ v i,2^ + C/t 3) r||^«||yi,2 (5T ) (3.25) 

i=l 

holds provided that u has the regularity demanded by the right-hand sides of the estimates 
and 5 € [0,1). IfTis an element with h\p = h-2,T < Ch 3) T an d St has zero distance to 
both a singular vertex and a singular edge then the local interpolation error estimate 

2 

\u — Dhu\ H l( T ) 

u \\v 1 ' 2 (s T ) ^1,T^3,t||^3w||^ 1 . 2 (5 t ) (3.26) 

i=l 

hold provided that u has the regularity demanded by the right-hand sides of the estimates 
and j3,6 € [0,1), + 6 <1. 

Proof. As in the proof of Lemma 3.7 we distinguish between the derivatives d 3 DhU and 
the derivatives along directions perpendicular to the X3-axis. From Lemma 3.3 we obtain 
by using the triangle inequality and 1 5*^ | — 1 ^ 2 1 1 c?3 f - 1 1 z, 1 ( S" r ) < ll^3 n llL 2 (T) 

\\d 3 (u-D h u)\\ L2{T) < \\d 3 u\\ L2iT) + C\ST\- 1/2 h a \\D a d 3 u\\ Ll{ST) 

<C\\d 3 u\\ L2{ST) + C\S T \- 1/2 h a \\D a d 3 u\\ L1{ST) . 

\a\=l 

For the estimate of diD^u, i = 1,2, we use Lemma 3.6, from which we conclude that 
\\di(u - D h u)\\ L 2( T ) 

h 3 ,T i 

\o-m\ L i { s T ) + 

These two estimates can be summarized by using hi t < Ch 3t T to 
\u - D h u\ H i(^t) 

<C\u\ hHSt) + C\S t \- 1/2 \^\\d 3 u\\ LHST) + £ h T \D a u\ wl , HST) ) . (3.27) 
It remains to estimate the terms against the weighted norms. Firstly, we have 

2 

\u\h HSt) < £ \\R l -^ l ~ s ■ R^e^diuW^T) + WR 1 -^ ■ Rf*- x e- x dzu\\v {T) 
i=i 

2 

<y^m&xR 1 ~ l3 1 ~ s \\diu\\ J ri,2 / r, x + m&xR 1 ^ 13 \\d 3 u\\,,i,2<r, v 
i=i 



< C\u\ HHSt) +C\S T \- 1 ' 2 [ ^ll^bi^) + £ ^Tl^kM(5 2 
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With Rt-PO 1 - 6 = r^R-PR 6 < r l - p - s R 5 < Ch{~ T ~ S h 3 3T (where we used the assump- 
tion f3 + 5 < 1) and R 1 ^ < h\~ T P < h~^,h 3 . T we derive 

2 

\u\ H 1( St ) < C h \~T~ &h i,T^2 W diU \\v±'1{S T ) + Ch l,T h 3,T\\ d 3 u \\ V ^(s T )- 

i=l 

With R 1 ~ s 8 1 ~ s = r 1 ^ 5 < r 1 ^ 6 < Ch\^ (using that the exponent is positive) we derive 
also 

2 

\u\ h i [St ) < Ch{~ T S Y^ WdMlv, 1 f (s T ) + Ch zAdM\ v ^{s T y 
1=1 

Secondly, for T intersecting the singular edge, but no singular vertices, we have 
^-\\d 3 u\\ L i {ST) < ^\\d 3 u\\ v i,2 {ST) \\r\\ L 2 iST) < /is.tI^tI 1 / 2 !!^!!^^). 

If T has also a singular vertex, then we have with RP~ 1 9~ 1 = R^r^ 1 

j^\\d3u\\ L i iST ) < ^\\d 3 u\\ v i : 2 (ST) \\R-Pr\\ L 2 {ST) < h 3yT h~ P T \S T \ l/2 \\d 3 u\\ v ia {ST) 
where we used that 

II^|| l2(St) < Wr 1 -^^ < Ch\J\S T \ 1 ' 2 (3.28) 

which can be obtained by integration. The second derivatives in estimate (3.27) are 
treated in a similar way. For i = 1, 2, 3 we get 

II^3'«||l1(S , t ) < IIl 2 (5 t ) 

For i, j = 1,2 and supposing that T does not have singular vertices we have 

h i)T \\d i3 u\\ Ll{ST) < h hT \\R- s 8~ s \\ LHST) \\R s e 5 d ijU \\ L 2 iST ) 
< h{^\S T \^\dM Vs y iST) , 

where we used again an argument as in (3.28). If T has a singular vertex, then 

KrWdijuWLUST) ^ hx, T \\R~ p 0~ \\ L 2 {St) \\R^ s d ljU \\ L 2 {ST) . 

But, R-P0- & = R-P +5 r- & < R 5 r-P- & < h^r-^ 5 , and so, since + S < 1, a similar 
argument as in (3.28) give us 

\\R-^- 5 \\ L 2 {St) < h{ T \\r-P- s \\ L 2 {ST) < Ch-^hySrl 1 / 2 . 

Hence we have 

h i,T\\diju\\ L i (ST) < Ch\^ s h^ T \S T \ ll2 \\diu\\ v i,2 {ST) . 
Therefore, the desired estimates are proved. □ 
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Theorem 3.10 (global interpolation error estimate). Let u be the solution of the 
boundary value problem (1.1) with f £ L 2 (fl), and let uj,ur be the functions obtained 
from the splitting (2.1). Assume that the refinement parameters m and i/£ satisfy the 
conditions 

M<\i £ \ (3.29) 
vi<>& + \, (3-30) 

v + h ( x ' l) -D >i > (3 - 3i) 

t = 1, . . . , L. Then the global interpolation error estimate 

\u R - D h u R \ HHKl) < Ch\\f\\ L 2 {Ae) (3.32) 

is satisfied. 

Proof. The proof can be carried out following the lines of the proof of Theorem 5.1 in 
[4] with the setting p = 2. Note that only a finite number (independent of h) of the St 
overlap at any point. □ 

Remark 3.11. The refinement conditions (3.29)-(3.31) were discussed in [4] already: 
The conditions (3.29) and (3.30) balance the edge and vertex singularities. The third 
condition, (3.31), follows from (3.30) in the case m = vi; only in the case fj,£ < V£ 
it imposes a condition between \i£ and vg limiting the anisotropy of the mesh. For the 
Fichera example treated in Section 4 we have Av ~ 0.454 and Ac = 3 . With the choice 
U£ = 0.9 the conditions (3.29) and (3.31) imply the choice 0.414 < /j,£ < |. For U£ = 0.8 
we would get the weaker condition 0.184 < /j,£ < |. 

Note also that in the absence of singularities we have set a| = o° and/or Ay = 00. 
In these cases we can set ni = \ and/ or V£ = 1. 

Corollary 3.12 (H 1 and L? finite element error estimate). Let u be the solution of 
the boundary value problem (1.1), and let Uh be the corresponding finite element solution 
on a finite element mesh as constructed in Section 2 with grading parameters satisfying 
the conditions (3.29)^(3.31). Then the discretization error can be estimated by 

\\u - Ufclljji(n) < C^ll/lli2(n): (3-33) 
\\u-u h \\ L 2 {n) <Ch 2 \\f\\ LHn) . (3.34) 

Proof. We choose = ui+DhUR in estimate (1.4) and observe that u—Vh = uh—D^ur. 
With Lemma 3.10 we obtain the estimate (3.33). The L 2 -error estimate can be derived 
by the standard Aubin-Nitsche method. □ 

Remark 3.13. A trivial conclusion from (3.33) is the stability estimate 

IKIIifi(fi) < C||/Ili2(n) (3-35) 

which we will need in Section 5. 
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Figure 4: Modification of macroelement of type 4 



Remark 3.14. In macroelements of type 4 with fi£ = vg < 1, Apel and Nicaise suggested 
in [4] the use of a more elegant refinement strategy as depicted in Figure 4- Our proof 
cannot be transfered to this kind of mesh immediately since there may be elements T 
where St is not prismatic as it was exploited in the proof of Lemmas 3.5 and 3.6. We 
conjecture that the assertion still holds but do not pursue this further in this paper. 

4. Numerical test 

As in [4] we consider the Poisson problem (1.1) in the "Fichera domain" Q, := (— 1, l) 3 \ 
[0, l] 3 and choose the right-hand side / = 1 + R~ 3 / 2 ln~ 1 (i?/4) which is in L 2 (0) but not 
in L p (£l) for p > 2. For this problem we have A„ ~ 0.45 for the concave vertex [29] and 
A e = ^ = | for the three concave edges. All other edges and vertices are non-singular. 

This boundary value problem was solved on quasi-uniform and on graded meshes 
with our refinement strategy using \i = v = 0.5 < min{A e ,A„ + where types 1, 
2 and 4 occur. Additionally we include the strategy where the macros of type 4 are 
replaced by type 5, compare Remark 3.14. Pictures of such meshes can be found in [4]. 
The refinement strategies and an a posteriori error estimator of residual type [31] were 
implemented into the finite element package MooNMD [19]. The estimated error norms 
are plotted against the number of unknowns in Figure 5. We see that the theoretical 
approximation order h 1 ~ iV -1 / 3 from Corollary 3.12 can be verified in the practical 
calculation for both refinement strategies. The error with the second strategy is slightly 
smaller. We denoted by iV the number of nodes. 
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Figure 5: Plot of the estimated error against the number of unknowns. The labels at the 
curve denote the estimated convergence order in terms of h ~ iV -1 / 3 . 

5. Discretization error estimates for a distributed optimal 
control problem 

Hinze introduced the variational discretization concept for linear-quadratic control con- 
strained optimal control problems in [17]. We follow here this concept in a special case. 
Consider the the optimal control problem 

rl \ 1 II ||2 II II 2 

(y,u)&H^(Q)xU !id I y ' Z ^ > 

where the state y G Hq(VL) is the weak solution of the Poisson problem 

— Ay = u in $7, y = on dCl, (5-1) 

and the control u is constrained by constant bounds a a ,u\, G R, this means that the set 
of admissible controls is defined by 

f/ ad := {u G L 2 (0) : u a < u < u b a.e. ft}. 

The regularization parameter a is a fixed positive number and y& G L 2 (Q) is the desired 
state. It is well known that this problem has a unique optimal solution (y,u). There 
is an optimal adjoint state p G Hq(Q), and the triplet (y,u,p) satisfies the first order 
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optimality conditions 

(Vy, Vv) L 2 {n) = (u,v) L 2 (n) 
(Vp, Vv) L 2 {n) = (y- y d , v) L 2 {n) 
(au + p,u- n) L 2( Q ) > 

With the variational discretization concept the approximate solution is obtained by 
replacing Hq(Q) by a finite element space Vh C Hq(Q) and searching {yhi^hiPh) £ 
V h x [/ ad x V h such that 

(yVh, V«ft)i2(n) = (Uh,Vh)L2(Q) V v h G Vft, 

(Vp h ,Vv h ) L 2^ = (y h - yd,Vh)L 2 (n) Vufc G Vh, 
(au h +ph,n - u h ) L '\n) > Vtt G [7 ad . 

Note that the control space is not discretized; nevertheless Uh can be obtained by the 
projection of —ph/a onto U ad , see [17]. The discretization error estimate 

II" - u h \\ L 2 {n) + \\y - yh\\L*(si) + \\P - Ph\\L*{si) < Ch2 (Il«lli2(n) + I|y<illi^(n)) 

can be concluded from (3.34) and (3.35), see [17, 7]. With the proof of Corollary 3.12 
we have established this result for anisotropic discretizations of the state equation (5.1) 
in the case of three-dimensional polyhedral domains. 



G H^(Q), 
Vv G H^(Q), 
Vn G U aA . 



6. Discrete compactness property for edge elements 

The Discrete Compactness Property is a useful tool to study the convergence of finite 
element discretizations of the Maxwell equations, both for eigenvalue and source prob- 
lems. It was first introduced by Kikuchi [20] and proved for Nedelec edge elements of 
lowest order on tetrahedral shape regular meshes. We refer to the monograph by Monk 
[23] and the references therein for further analysis on isotropic meshes. The property 
was also analyzed on anisotropically refined tetrahedral meshes on polyhedra for edge 
elements of lowest order by Nicaise [24] (excluding corner singularities) and by Buffa, 
Costabel, and Dauge [13]. 

Lombardi [21] extended this result to edge elements of arbitrary order, also including 
corners and edge singularities. The proof is based on two tools: 1) interpolation error 
estimates for edge elements on meshes satisfying the maximum angle condition, and 
2) interpolation error estimates for a piecewise linear interpolation operator defined on 
W 2 ' p (Q)nHQ (O), p > 2, preserving boundary conditions. For the latter, the Lagrange in- 
terpolation was used (implying p > 2) together the results of Apel and Nicaise [4], giving 
some artificial restrictions on the grading parameters defining the allowed anisotropically 
graded meshes. Using now estimate (3.33) of Corollary 3.12 we can extend the result of 
[21] allowing little more general meshes. 
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In what follows we define a family of edge element spaces and introduce the DCP for 
this family. We refer to [21] for further definitions and notation. First we introduce the 
divergence-free space 

I = {v£ J3o(curl, 0) : divv = on ft} . 

Then we introduce discretizations of this space where the divergence-free condition is 
weakly imposed. Let I be a denumerable set of positive real numbers having as the 
only limit point. From now till the end of this section, we assume that h El. For each 
h, let Tii be the mesh on the polyhedron ft constructed in Section 2. Given an integer 
k > 1, let Xh be the space defined as 

X h = {v h , G #o(curl,n) : v h \ T E A/" fc (T)VT G %, (Vp hj v A ) = 0Vp h G S h } 

where jVfe(T) is the space of edge elements of order k on T, and 

S h = { Ph G Hfcn) : p h \ T G V k (T)VT G T h } . 

We say that the family of spaces {Xh}hei satisfies the discrete compactness property if 
for each sequence {v^j^gj, JCI, verifying for a constant C 

v h G Xh, V/i G J, 
ll v h||i? (ctirl,Q) < C, V/t G J, 
there exists a function v G X and a subsequence {v/j n } ne N such that (for n — > oo) 

v hn ^ v in L 2 (Q) 

v h„ ~^ v weakly in Hq(cut1, f2). 

Theorem 6.1. If the grading parameters defining the meshes Th satisfy the conditions 
(3.29)-(3.31), then the family of spaces {X^^o verifies the discrete compactness prop- 
erty. 

Proof. Follow exactly the arguments used to prove Theorem 5.2 of [21] taking into ac- 
count that the inequality (4.21) of that paper is now a consequence of estimate (3.33). □ 

A. Proof of trace inequalities 

Lemma A.l. Let P be a triangular prism with vertices Vi, i = 1, . . . , 6, where the face 
V1V2V3 is opposite to the face v^v^vq, and where the edges V1V4, V2V5, and v^vq are parallel 
to the x^-axis, see Figure 6. Denote by F the face V1V2V3. Then for all v G W 1,P {P), 
p G [1, 00), we have 

\H p lp {f) < ^ • ^3 1 (JMIlp(p) + K\\ d M\ p LP{ P)) , 

where /13 is length of the shortest vertical edge, and^f is the angle between the x\X2-plane 
and the plane containing the face F. The constant C reg depends only on the minimum 
angle of the face F. 
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Figure 6: Illustration of the prism 



Proof. We can assume v% = (0,0,0) and V4 = (0, 0, /13). Suppose V2 = (0,2,02,02), 
v 3 = (03,03, £3)- Let s,t such that 

a 2 s + b 2 t = c 2 
a 3 s + b 3 t = c 3 . 

It is clear that there exist such s and t since vi, V2, and v 3 do not lay on one line. Then 
the map f(x) = Bx, with 

10 
D = I 1 

s t 1 

sends P to P where P is a prism with three vertical edges and some of its vertices are 
v\ = (0, 0, 0), V2 = (02, 02, 0), V3 = (03, 63, 0) and £4 = (0, 0, /i 3 ). Let F be the face v\V2V 3 
of P. 

Let v be defined by ■y(x) = v (x) ii x = Bx. Then we have 

Now, if Q is the right prism with vertices Vi, . . . , V4, (02, b2, h 3 ) and (a 3 , v 3 , /i 3 ), then we 
have using a trace inequality on Q and noting that Q C P that 



'LP(F) 



'LP(Q) 



W\V.-^ + hl\\d,vf Tvl ~ 



'Lp(P) 



(Q) 
Lp(P) 



with C p depending only on p. Therefore, we have 



a 



iLP(F) 



COS 7 



^l'< p(i5) + 



c, 



cos 7 



V W P LP(P) + h 3\\ d 3v\\ P LP{P) 



where we used that d 3 v(x) = d 3 v(x). Since \B\ = 1 we obtain the desired result. 



□ 
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Figure 7: Illustration of the notation used in Lemma A. 2. The dotted lines indicate the 
prism Q n , dashed lines the parallelogram P n while the triangle £(z) is hatched. 
Note that a(z) = n F n . 

Lemma A. 2. Let T be an anisotropic element with the node n on the singular edge and 
let o~ n be a short edge. Let P n C St be a parallelogram of maximal area having a n as 
an edge and another edge on the singular edge, see Figure 7. And let F n the face of St 
containing P n . Then \P n \ > C|-F n |, and for all v 6 W 1,1 (St) we have 

|MlLi(P n ) < C\F n \\S T \~ l {\\v\\ L i {St ) + l«l,T|||fla lir «|Ui(S T ) + \s2,T\\\d S2tT v\\ L l( ST) ) . 

where si ; t and S2,r are two short edges ofT. 

Proof. The inequality \P n \ > C|-F n | follows from our assumptions on the mesh, in par- 
ticular from the comparable length of opposite edges of F n . For proving the estimate 
choose the coordinate system such that n = (0, 0, 0). 
Assume first v is regular. We have 

||«IUi(Pn) ^ C / \v({0,0,z) + ta n )\ dtdz 

Jo Jo 

/ \v\ ds dz 

Ja(z) 

where cr(z) is the segment parallel to a n and with the same length and passing through 
(0, 0, z). If is the triangle contained in St having a(z) as an edge and being parallel 
to the bottom face of St, then since we can assume v\^/ z \ is regular (because v is it), by 
a trace inequality we have 

/ M < [ (\v\ + \s 1 , T \\d SlT v\ + \s 2 ,T\\d S2T v\) 

Ja(z) I? I J£(z) 
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where |si,t| and |s2,t| are the lengths of two small edges of T and |£| = |£(0)|. So we 
have 



\v\W{p n ) < C J7T / (M + \ s i,T\\ds hT v\ + \s 2 ,T\\d S2T v\] 

ISI Jo J£(z) 

\F n \ f hi ' p « f 
- C T^1 / (\ v \ + \ s hT\\d SltT v\ + \s 2 ,T\\d S2iT v\ 

\ b T\ Jo j£(z) 



v\) 



- C W\ I (\ v \ + \ s hT\\d Sl , T v\ + \s 2 ,T\\d S2 . T 

\ b T\ jQ n 

- C \^\ I (1^1 + \ s 1,t\\ 9 si,t v \ + \ s 2,T\\d S2:T v\) 

\ b T\ Js T 

where Q n is the prism formed by the union of £(z) with z G [0, ^3,p n ] that is contained 
in St- 

live W 1,1 (St), let {-Ufcjfc be a sequence of (7°° functions converging to v in W 1,1 (St)- 
For each /c we have 

||«fc|Ui(P„) < ClFnWSr^ 1 (||«fc|Ui(s T ) + l*i,r|||3 aiiT u fc || £ i (ST ) + |« 2 ,t| \\d S2tT v k \\ L i { s T )) ■ 
Now, the proof concludes by taking limit k — > oo. □ 
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